Three- Wave Modulational Stability and Dark 
Solitons in a Quadratic Nonlinear Waveguide with 
Grating 

Arthur Gubeskys and Boris A. Malomed 

Department of Interdisciplinary Studies, School of Electrical Engineering, Faculty of 
Engineering, Tel Aviv University, Tel Aviv 69978, Israel 

Abstract. We consider continuous-wave (CW) states and dark solitons (DSs) in a 
system of two fundamental-frequency (FF) and one second-harmonic (SH) waves in 
a planar waveguide with the quadratic [x^ 2 -*] nonlincarity, the FF components being 
linearly coupled by resonant reflections on the Bragg grating (the same model is known 
to support a great variety of bright solitons). We demonstrate that, in contrast with 
the usual situation in spatial-domain \^ models, CW states with the phase shift 7r/2 
between the FF and SH components are modulationally stable in a broad parameter 
region in this system, provided that the CW wavenumber does not belong to the 
system's spectral gap. Stationary fundamental DSs are found numerically, and are also 
constructed by means of a specially devised analytical approximation. Bound states 
of two and three DSs are found too. The fundamental DSs and two-solitons bound 
states are stable in all the cases when the CW background is stable, which is shown by 
dint of calculation of the corresponding eigenvalues, and verified in direct simulations. 
Tilted DSs are found too. They attain a maximum contrast at a finite value of the tilt, 
that does not depend on the x^ 2 ^ phase mismatch. At a maximum value of the tilt, 
which grows with the mismatch, the DS merges into the CW background. Interactions 
between the tilted solitons are shown to be completely elastic. 
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1. Introduction 

Optical dispersive and diffractive media with quadratic [x^] nonlinearity are well known 
for their potential to support various types of solitons, which have been a subject of 
intensive studies [1,2]. In most cases, x^ solitons are experimentally observed in the 
spatial domain, as the small size of available samples makes it difficult to accommodate 
the dispersion length of temporal solitons. Nevertheless, using special techniques, such 
as tilted wave fronts, it is possible to induce a strong artificial dispersion and thus create 
temporal solitons in available x^ optical crystals [3]. 

Another challenge for the experiment was creation of x^ dark solitons (DSs) 
and their two-dimensional counterparts (optical vortices). For the first time, DSs 
in x^ models were considered in the cascading limit, which corresponds to a large 
phase mismatch between the fundamental-frequency (FF) and second-harmonic (SH) 
waves [4]. In this limit, the x^ model reduces to the nonlinear Schrodinger (NLS) 
equation with the Kerr [x^] nonlinearity. Beyond the cascading limit, a single analytical 
solution [5] , and a continuous family of numerical solutions for fundamental and twin- 
hole DSs [6, 7] were found. All these solitons are unstable in the spatial domain, due to 
the modulational instability (MI) of the CW (continuous-wave) plane-wave background 
which supports the DS. The only possibility to avoid the MI was found in the temporal 
domain, in the case when the group- velocity-dispersion coefficients have opposite signs at 
the FF and SH [8-10] (in fact, such a case is quite realistic in the temporal domain [11]). 
Nevertheless, effectively stable spatial-domain vortices, supported by a background of 
finite extension, were experimentally created in the case of negative phase mismatch 
between the FF and SH waves [12]. In the latter case, the MI had the character of a 
convective instability, so that perturbations were carried away from the vortex's core 
faster than they grew due to the MI. 

Additional possibilities for creation of solitons are offered by a combination of x^ 
nonlinearity with Bragg gratings (BGs). The grating may be introduced in both the 
time-domain and spatial-domain models, but in the latter one it is a more realistic 
feature, amounting to a system of parallel "scratches" drawn on the surface of a 
quadratically nonlinear planar waveguide. In that case, the generic x^ model is a 
four-wave one, as both the FF and SH components include two counter-propagating 
waves (detailed descriptions of the model can be found in the reviews [1] and [2]). The 
MI of CW solutions was studied in the four-wave model, with a conclusion that they 
are unstable in most cases. Stability regions for the CW states were also found, close to 
the transition to the NLS limit, but only in a rather exotic case when the BG strength 
at the SH is larger than at the FF [13]. A stable dark bi-soliton in this model was found 
numerically near the band edge [14]. 

In this paper, our aim is to investigate the MI of CW states, and, in the case 
when they are stable, to find DSs in a three-wave spatial-domain system, in which 
two FF components are linearly coupled through the resonant scattering on the BG. 
Simultaneously, the quadratic nonlinearity couples the two FF components to the SH 
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wave. The propagation direction z is parallel to the "scratches" which form the BG. 
Then, the above-mentioned couplings can be achieved by choosing FF carrier wave 
vectors to be of equal magnitude and making opposite angles with the z axis, while 
the SH wave vector is parallel to it. Thus, only the FF components are affected by 
the grating, while the SH wave is subject to diffraction. This model was introduced in 
Ref. [15] in the context of bright solitons. Stable gap solitons, as well as bright embedded 
solitons (which exist inside the continuous part of the system's linear spectrum [16]) 
have been found in this three-wave model. 

The paper is structured as follows. The model is formulated in sectional where we 
also give an estimate for relevant values of the physical parameters. In section El we 
identify two types of CW solutions and analyze their modulational stability. It is shown 
that, on the contrary to other x^ models in the spatial domain, in the present case one 
type of the CW solutions is stable in a broad range of parameters. DSs supported by the 
stable CW background are considered in section 0] For stationary fundamental solitons, 
we find an approximate analytical and direct numerical solutions. Tilted solitons (with 
a slant relative to the propagation axis z), as well as bound states of two DSs, are found 
too. Section presents results for the stability and interactions of the DSs. It is shown 
that both the fundamental solitons and their bound states are stable in the whole region 
where the respective CW is modulationally stable. The work is concluded by section |3 

2. The model 

We consider a quadratically nonlinear planar waveguide with a spatial BG in the form 
of scores directed parallel to the propagation direction z. As shown in Ref. [15], spatial 
evolution of two components U\ and u 2 of the complex FF field, and the complex SH 
field W3 obey the following system of normalized equations: 



Here z and x are, respectively, the propagation and transverse coordinates, the subscripts 
stand for partial derivatives, the asterisk means for the complex conjugation, q is the 
phase-mismatch parameter, and D is an effective diffraction coefficient, defined with 
regard to the fact that Bragg- reflectivity and x constants are normalized to be 1. We 
notice that the complex conjugation applied to Eqs. (JH), and the change of the notation 
u\ — > Ui,u 2 — > —u 2 ,u^ — > u 3 , q — > — q transforms D to —D, therefore we may assume, 
without loss of generality, that D is always positive (while q may be positive, negative, 
or zero). 

Stationary solitary-wave solutions are looked for as 



% (u 1 ) z + i {u 1 ) x + u 2 + u 3 u* 2 = 0, 
i (u 2 ) z - i {u 2 ) x + Mi + u 3 u\ = 0, 
2i (u 3 ) z + D (u 3 ) xx + u x u 2 - qu 3 = 0. 



(1) 



u i,2 = Ui )2 (x — cz) exp(ikz), u 3 = U 3 (x — cz) exp(2z/cz), 



(2) 
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where real k and c are the propagation constant and tilt of the corresponding optical 
beam in the (x, z) plane, and complex functions t/1,2,3 are to be found from the equations 

i(l - c)U{ - kUt + U 2 + U 3 U^ = 0, 

+ c)U' 2 - kU 2 + U x + U 3 U* = 0, (3) 

DU' 3 ' - 2icU' 3 - aU 3 + UxUi = 0, 
with a = 4k + q and the prime standing for d/d(x — cz). For symmetric (untilted) 
solutions with c = 0, one may set U\{x) = —U£(x) = U(x), U 3 (x) = V(x), where the 
function V(x) is real, which yields a simplified system, 

iU' - kU - U* - VU = 0, 

DV" -aV - \U\ 2 = 0. ( ' 

In physical units, the same value of the power density of the CW background in 
the LiNb03 planar waveguide that was used in the experiments with bright spatial 
solitons [17], i.e., ~ 20 W//zm at the pump wavelength ~ 1.3 /im, may be assumed, 
in combination with the BG reflectivity ~ 1 cm -1 , which is a typical value for weak 
gratings. Then, the physical value of the wavenumber detuning corresponding to k > 1, 
i.e., just above the upper edge of the spectral gap, which is shown below to be appropriate 
for DSs, is about ~ 1 cm -1 . The range of the parameters a and D relevant to the 
experiment can be determined as in Ref. [15], where the normalized parameters of the 
model were related to physical ones. In particular, a may take values in a wide range, 
4 < a < 100. For example, the large phase mismatch used for the experimental generation 
of vortices in a media in Ref. [12] corresponds to a ~ 60. Realistic values of the 
effective SH diffraction parameter are 0.1^ D< 1, depending on the angle between the 
FF wave vectors and the z axis. The angle must be, however, small enough for the 
applicability of the paraxial approximation, which is implied in Eqs. (JIJ. 



3. Modulational stability of CW solutions 

Equations (0J) give rise to two types of CW solutions. The first of them is a real one, 
with in-phase FF and SH components, 

U = U = ±y/a(k + l), V = V = —(k + 1), (5) 

which exists provided that a(k + l) > 0. Below, it will be referred to as a PR (pure-real) 
solution. In the other solution, the phase of the FF component is shifted by 7r/2 against 
the SH, 

U = U = ±iyja(k-l), V = V = l-k. (6) 
This solution, with a purely imaginary FF field, will be accordingly called a PI one (of 
course, the SH is real in it). The PI solution exists provided that a{k — 1) > 0. 

To investigate the modulational stability of the solutions, we consider a perturbed 
one in the form of 

u\ = [u[°^ + af'* exp(— iXz + iQx) + bi exp(i\*z — iQx)] exp(ikz), 

u 2 = [U^ + CL2 exp(— iXz + iQx) + b^ exp(i\*z — iQx)} exp(ikz), (7) 

u 3 = [U 3 QS) + a 3 exp(—iXz + iQx) + b^ exp(iX*z — iQx)] exp(2ikz), 
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where u[^,3 stand for the amplitudes of the CW solutions defined above, 0^2,3 > &i°2,3 an d 
Q are infinitesimal amplitudes and arbitrary real wavenumber of the perturbation, and A 
is the eigenvalue to be found. Substitution of these expressions into linearized equations 
leads to a resolvability condition for jj, = A 2 , in the form of a cubic equation 

^ + T 2 (Q)fi 2 + T 1 (Q)fi + T (n) = 0. (8) 

Obviously, the stability requires that \x must be real and positive for all Q. 

We start the analysis for the case of D = 1. Then, the coefficients in Eq. (jHJ are 

T 2 = -ifi 4 - (la + 2)Vl 2 - 2ak - \a 2 + 4k - 2a, 

T 1 = f ft 6 + (a + 1 - k)tt 4 + (\a 2 -2ak + a- ak 2 - Ak 2 - Ak)VL 2 - ka 2 - 8ka - 8k 2 a, 
T = -\n 2 (-Aka + att 2 - 4a + Q 4 ){-Ak 2 - Ak + a + Q 2 ), 

for the CW solution of the PR type. For the Pi-type solution, they are 

T 2 = -ifi 4 - (fa + 2)n 2 - 2ak - \a 2 - Ak + 2a, 

T 1 = ifi 6 + (a + 1 + A;)fi 4 + (±a 2 + 6aA; + 3a - aP - Ak 2 + Ak)n 2 + ka 2 - 8ka + 8k 2 a 
T = \n 2 (Aka + ail 2 - 4a + VL A ){Ak 2 -Ak-a-tt 2 ). 

Equation (jHJ) was solved numerically (an analytical solution is formally available, 
but it takes an intractably complex form). The result, in the form of the stability region 
of the PI solution, is displayed in Fig. ^ Notice that the solution is always unstable if 
its wavenumber belongs to the gap of the FF subsystem, which is k 2 < 1. 

The CW solutions of the PR-type solutions are always modulationally unstable. 
However, for large values of |a| (actually, in the case of large phase mismatch), the 
instability band of the perturbation wavenumber Q becomes very narrow. For example, 
it is 8.69 < < 8.71 for a = —100 and k = —1.2. Therefore, the instability generated 
by a random perturbation will grow very slowly in that case, and, for practical purposes, 
the CW state may be regarded stable. 

The above results pertain to D = 1 in Eqs. With other values of D, we have 
found that, for the PI- type CWs, the increase of D results in shrinkage of the stability 
area, and the decrease of D results in its expansion, see the dotted and dashed curves 
in Fig. ^ The PR-type CW solution remains unstable at any D. 

4. Dark solitons 

4-1. Stationary solitons 

The availability of the modulationally-stable CW background suggests a possibility of 
existence of stable DSs. As concerns search for stationary symmetric (c = 0) solutions, 
it is obvious that Eqs. (J3J are invariant with respect to the transformation U — » U / VD, 
V — > V, a —>■ a/D, and D —> 1, which makes it possible to set D — 1, varying only a 



Modulational Stability and Dark Solitons in Nonlin. Waveguide with Grating 6 



in the stationary equations. Further, we substitute U(x) = p(x) + ig(x), which yields 
equations for the real functions p, g and V, 

j/ = g(k-l+V), 

g' = - P (k + l + V), (11) 
V" =p 2 + g 2 + aV. 

The SH component can be sought for as V — Vq + V, where Vo is the CW solution 
(JHJ) of the Pi-type. Then, for the solutions supported by the Pi-type CW background, 
we obtain, from Eqs. (fTTj) . 

p' =gV, g' = -p(2 + V), ^ (12) 
V " = p 2 + g 2 + a(l - k) + aV. (13) 

Assuming, for the time being, V < 1, we approximate the second equation in (fT2"J) by 
g' = —2p (the validity of this assumption will be considered later). Then, we look for a 
fundamental DS solution in the form 

g = \U Q \ tanh(/fc), p = p sech 2 (/3x), 

V = Vb + y Q sech(/3x), 1 ' 

where U and V are the same as in Eqs. © (hence the solution is matched to the 
CW background at x — > ±oo), while the constants po, Vq and j3 are to be found. Note 
that, according to Fig. [IJ Vo = 1 — k is negative in the region of stability of the CW 
background, hence the expression (JT3jl with Vo > assumes that |V(x)| has a minimum 
at x = 0. Substitution of the ansatz (JHJ) into Eqs. (|T2"j) and (fT3|) shows that it solves 
the first two equations, but not the last one. To identify an approximate solution, we 
then demand the latter equations to be satisfied only at the central point, x — 0. This 
approach leads to a system of algebraic equations, which can be readily solved to yield 
the following results: 



Vj = l ( 1 _ 14/ 4)U--^^ i -l), (15) 



(3=Vl / o, Po = -dlWVo (16) 



provided that Vo > and 1 — (8Vo/a) (1 — Vo/4)~ 2 > 0. It is straightforward to check 
that the latter conditions are satisfied in the region of interest (a > 0, k > 1), where the 
CW solutions of the Pi-type are modulationally stable, see Fig. [TJ 

The above approximation assumed Vo < 1. To examine the validity of this 
assumption, we note that the expression (|T5|). considered as a function of a, attains 
the maximum value (Vo) = 4Vo/ (Vo — 4) = A(k — l)/(3 + k) at a = oo, hence the 

V / max 

assumption is formally justified for small (k — 1) (near the upper edge of the spectral 
gap) if a (i.e., the phase mismatch) is large, and for larger (k — 1) if a is smaller. 

According to Eqs. (JHJ), the parameters po and Vo determine the inverse DS 
contrasts <J u ,v, which may be defined as ratios of the amplitudes at the soliton's center 
to those in the CW background: a u = |po/^o| an d cr„ = |Vo/Vo| for the FF and SH, 
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respectively. These parameters belong to the interval < a UjV < 1, the limiting cases 
of a u>v = and a U)V = 1 corresponding, respectively, to the black soliton, and to the 
unmodulated CW solution. In particular, Eqs. (fTK|) and (fi7)|) predict the black soliton 
at a = 0. 

The dependences of the inverse contrast on a, as obtained from the above 
approximation and from a direct numerical solution, are shown in Fig. |21 (two cases 
shown in this figure adequately represent the general situation). The results demonstrate 
that the analytical approximation is, generally, quite accurate, although numerically 
found solitons are never completely black. With the increase of k, the FF component 
becomes "grayer" , while the SH one gets darker. 

Typical examples of the direct comparison between the DS shape, as predicted by 
the analytical approximation, and as generated by the numerical solution, are displayed 
in Fig. El For k = 1.2 and a = 20, the approximation is so close to the exact solution 
that they are indistinguishable. For k = 2.4 and a = 40, the accuracy deteriorates, but 
the agreement is still fairly good. 

Further numerical search for DS solutions supported by the modulationally stable 
CW background has revealed the existence of bound states of the DSs, whose shape is 
similar to that of unstable bound states which were reported Refs. [7] and [14] (however, 
the bound states may be stable in the present model, see below). Typical examples of 
twin-hole and tri-hole DS complexes are shown in Figs. lUand^l 

4-2. Tilted solitons 

The general ansatz (j2J) for steady-shape solutions admits tilted DSs with c ^ 0. We 
have found such solitons by numerical search in the (k, a) plane. The above-mentioned 
scaling, which allowed us to impose the normalization D = 1 for the c = solitons, 
does not hold for the tilted ones. Nevertheless, we set D — 1 in this case too, as 
the additional problem of scanning the parameter space augmented by the effective 
diffraction coefficient D is too complex. 

For the tilted DSs, the amplitudes of the two FF components Ui and u<i are no longer 
equal, and the SH field is not real; accordingly, one should use Eqs. (jHJ) instead of Eqs. 
(J3J). Numerical solution for the tilted ("moving") solitons gives rise to the dependence of 
the inverse contrasts in the three components ( "forward" and "backward" FF ones, and 
the SH component), defined as above, on the slope (tilt) c. A typical result is shown 
in Fig. |H1 For the forward (backward) slope of the DS, the contrast of the forward 
(backward) FF component initially increases and attains a maximum at some tilt c max . 
For |c| > |c max |, the contrasts of both FF components decrease until the DS reaches 
a point c = Cbif where o becomes equal to 1, i.e., the DS disappears, bifurcating back 
into the CW solution. For the forward-tilted soliton, the forward FF component of the 
soliton is always darker than its backward counterpart. 

The dependence of parameters of the tilted DS on a, i.e., on the phase mismatch, 
is shown in Fig. As is seen, the slope c max , at which the maximum contrast of FF 
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occurs, is virtually independent of a. However, the slope Cbif at which the tilted DS 
bifurcates back into the CW, i.e., the maximum slope past which the DS does not exist, 
increases with a. 

The dependence of the same parameters on k, 1.6., cLS cL matter of fact, on the 
detuning from the gap, is displayed in Fig. El As well as it was with the a-dependence, 
the largest FF contrast is independent of k, but the corresponding c max grows with k. 
The bifurcation value Cbif of the tilt is also a growing function of k. It was observed too 
that the maximum values of the contrast remain themselves constant with the variation 
of both a and k. 

As k becomes still larger than in the cases like that shown in Fig. |Hl viz., at k > 2.1, 
the behavior of the DS solution changes drastically. In this region, the numerical scheme 
ceases to converge close to the bifurcation point, therefore the bifurcation (merger of 
the DS back into the CW) cannot be identified. A typical example of the dependence 
of the DS contrast on the tilt for this situation is shown in Fig. |§] (cf. Fig. |)J). In this 
case, the numerical method does not converge for c > 0.22. 

The normalized slope c relates to its physical counterpart c p h ys , which is the ratio 
of the real coordinates, as c p h ys = pc, where p is the angle between the FF wave vector 
and the axis z. Since the p should be small for Eqs. (0) to remain valid, the physical 
slopes expected to be observed in the experiment are on the order of 10~ 3 — 10~ 2 . 

5. Stability and collisions of dark solitons 

The stability of the CW background is only a necessary, but not sufficient condition for 
the stability of DSs. We investigated their full stability by direct simulations, and also 
through a numerical solution of the eigenvalue problem for small perturbations. In the 
latter case, the perturbed solution was taken as 



Ui 
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exp(iAz) 


exp(ikz), 
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x\ k) + e 3 (x) 


exp(zAz) 


exp(2ikz), 





where u$>(x; k) represent an unperturbed DS solution, e n are eigenmodes of infinitesimal 
perturbations, and A is the corresponding eigenvalue, the stability requiring all A to be 
real. The expressions (fr7|) were substituted into Eqs. (0), and the resulting equations 
were linearized in e n . 

The problem significantly simplifies for the symmetric solitons with c = 0, when 
itj = — (^2°^ , and ?4 is real. In that case, we were able to find the full spectrum 
of the eigenvalues for both the fundamental DSs and their twin-hole bound states (see 
Fig. E}, on a sufficiently dense grid of points in (a,k) plane. As in the previous sections, 
the analysis was restricted to the case of D — 1, as scanning through the full three- 
dimensional parametric space including D is a very complex technical problem. As a 
result, it has been concluded that both the fundamental and twin-hole DSs are always 
stable, provided that their CW background of the PI type is modulationally stable. 
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The stability of the symmetric DSs was also verified in direct simulations of Eqs. (JTJ. 
The verification has shown that they are stable indeed, as predicted by the calculation 
of the perturbation eigenvalues. 

Computation of the stability eigenvalues for the tilted fundamental DSs is too 
hard to run it in a systematic form. However, their stability can be easily tested in 
direct simulations. The simulations reveal that, as well as their symmetric (untilted) 
counterparts, the tilted solitons are always stable if the CW background is stable. 

The stability of the DSs with c ^ suggests a possibility to consider "collisions" 
between them (in fact, intersection between two spatial solitons tilted in opposite 
directions). The collisions were simulated too, with a conclusion that they always 
look completely elastic: the DSs re-emerge unaffected after the collision, apart from 
a displacement of their centers, see an example in Fig. The character of the 

interaction remains the same for all values of a and k (at least, for the slopes \c\ > 0.1). 
The collision-induced shift of the centers is (quite naturally) slightly larger for smaller 
slopes. A typical X-j unction formed by the intersection of two dark solitons is shown in 
Fig. UU 

The collisions are, actually, quite similar to those observed in the integrable NLS 
model describing dark solitons in a defocusing Kerr medium. As the solitons collide, 
they repel each other, which is obvious from the contour-plot rendition in Fig. ^2 If 
is relevant to note that interactions between DSs in the integrable self-defocusing NLS 
equation are repulsive too (see, e.g., Ref. [18]). 

We also simulated a different type of interactions between the DS, when two solitons 
with c = are originally placed at some distance; obviously, this situation may be 
relevant to the experiment, and to possible applications. It has been observed that the 
solitons repel each other in this well. With very small radiation loss, they start to 

"move" (i.e., acquire a finite slope), as is shown in Fig. The eventual slope observed 
in this set of the simulations virtually does not depend on k, and slightly increases with 
a. 

6. Conclusions 

In this work, we have considered uniform CW states and dark soliton (DSs) in the 
system of three waves in a planar waveguide [with two fundamental-frequency (FF) and 
one second-harmonic (SH) components] coupled by the x^ interaction. The FF waves 
are also linearly coupled by reflections on the Bragg grating. This model was earlier 
shown to generate a great variety of bright solitons. 

We have demonstrated that, on the contrary to the usual situation in spatial-domain 
•^(2) models, in the present case CW states with the phase shift n/2 between the FF 
and SH waves are modulationally stable in a broad parameter region, provided that 
the CW wavenumber lies outside the system's spectral gap. Stationary fundamental 
DSs supported by the stable CWs were found numerically, and also by means of an ad 
hoc analytical approximation, which produced results that agree well with numerical 



Modulational Stability and Dark Solitons in x^ Nonlin. Waveguide with Grating 10 

findings. Bound states of two and three DSs were found too. The fundamental DSs, as 
well as dark bi-solitons, are stable in all the cases when the CW background is stable; 
this was demonstrated by means of calculation of the corresponding stability eigenvalues, 
and verified in direct simulations. 

We have also studied tilted fundamental DSs. It was found that they attain a 
maximum contrast at a finite value of the slope, that virtually does not depend on 
the phase mismatch; at a maximum value of the slope (which grows with the 
mismatch), the DS merges into a CW state. Interactions between tilted solitons were 
systematically simulated too, with a conclusion that the interactions are repulsive and 
completely elastic. 
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Figure 1: In the shaded region in the plane of the parameters (a = Ak + q, k), the CW 
solution (jUJ) of the PI type is stable for D — 1 [recall that k is the wavenumber of the 
solution, see Eqs. (J2J), and q is the phase- mismatch coefficient in Eqs. (0)]. The stability 
borders of the CW solution for D = 0.5 and D = 10 are shown, respectively, by the 
dashed and dotted curves. 
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Figure 2: The inverse contrasts a u and a v of the fundamental dark soliton supported by 
the modulationally stable CW background of the Pi-type. The solid and dashed-dotted 
lines show, respectively, the analytical prediction for a v and <r u ; curves of the same types 
connecting circles show numerical results for a v and a u . The panels (a) and (b) pertain 
to k — 1.2 and k = 2.4, and in both cases D — 1. 
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(a) k = 1.2, a = 20 (b) k = 2.4, a = 40 

Figure 3: Typical examples of the fundamental dark soliton. The solid and dotted 
lines show, respectively, the numerical solution and the analytical approximation: (a) 
k = 1.2, a = 20; (b) k = 2.4, a = 40. 
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Figure 4: A bound state of two dark solitons, found numerically at k — 2.8 and a = 40. 
In this figure and below, the labels FH and SH pertain, respectively, to the fundamental 
and second harmonics. 
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Figure 5: A typical example of a bound state of three dark solitons, found for k = 2.4 
and a = 20. 
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Figure 6: The inverse contrast of three components in the tilted ("moving") fundamental 
dark soliton vs. its tilt c for k = 1.2 and a = 20. At points c = Cbif, where all the 
three contrasts becomes equal to 1, the dark soliton disappears, merging into the CW 
background. 




Figure 7: The slope c m i n of the tilted dark soliton, at which the maximum contrast of 
the soliton's FF component is achieved (circles), and the slope Cbif, at which the dark 
soliton merges into the CW background (squares), versus the effective mismatch a, for 
k = 1.2. 
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Figure 9: The same as in Fig. |H1 but for k = 2.2. 



Modulational Stability and Dark Solitons in Nonlin. Waveguide with Grating 16 




Figure 11: Contour plot of two dark solitons collision showing the displacement of their 
centers. The slopes are c = ±0.1, while k = 1.2 and a = 40. 
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Figure 12: A typical example of the interaction between two dark solitons with the c = 
placed initially at some distance from each other. In this case, k = 1.2 and a = 20. 



